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see,  fair  example,  Thom,  Allen  Southvell,^  anA  Jcneou.^)  figurations  are  first  established,  and  then  contour  patterns  €U*c 


tlals  of  the  calcuOatlonal  method.  Third,  dlBousslons  are  given  of  the 


If  the  fluid  Is  Incompresslhle,  then 
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(2-15) 


Note  that  if 


check  whether  any  Instabilities  exist.  The  experiments  do  Indeed  Indi-  Growth  of  the  perturbation  was,  however,  observed  for  p  >  0.8,  Some 

cate  that  instabilities  will  result  if  the  condition  typical  examples  of  the  form  of  the  growth  are  displayed  in  Fig.  2.  Here 


permits  convergence  of  the  successive  ^proximation  solution  to  occur  in 
almost  half  the  time,  because  the  new  solution  is  displaced  less  from 


with  coinputer  methods,  is  to  tisc  any  available  advanced  values  of  nei^>- 


wlth  (5-1)  the  correspondence  with  the  continuity  equation  is  direct 
Z  of  the  kth  iterate  at  neighboring  points  terms  of  the  velocity  components,  whereas  the  interpretation  of  the 


(3-1)  that  are  uaed  may  differ  on  the  left  and  right  hoimdarlea,  result- 
ing  in  a  sli^t  difference  in  the  final  boundary  values.  If  a  suffidartly 
stringent  convergence  criterion  is  used,  these  differences  in  the  ri^t 


note  that  If  we  sum  the  vortlclties  along  the  positive  y  axis 


substitute  u  for  the  average  of  the  velocities  as  prescribed  by  (2-7). 


Again  note  that  averages  are  used  as  obtained  by  (2-7)  through  (2-10)  where  i^  is  the  x  index  for  the  ri^t  boundary, 

with  permutations  of  the  indices  as  required.  In  the  averages  the  re¬ 
spective  wall  velocities  are  subsituted  for  u,  i  and  u 


time  on  u)  using  (2-6)  everyvhere  except  at  the  obstacle.  '</e  then  let 
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separately  with  regard  to  the  aide  from  which  the  obstacle  Is  observed. 

Further,  tests  would  have  to  be  ^Ued  In  evaluating  these  vortlcltles,  current  mesh  used  are  such,  however,  that  an  Inflow  <,f 

using  (2-4),  such  that  proper  velocity  values  would  be  used.  Velocity  Interpreted  to  represent  a  Jet,  since  obstacles  can- 

values  on  opposite  sides  of  the  obstacle  certainly  cannot  be  used  for  placed  sufficiently  far  downstream. 


prescription  affects  the  ca^.culation 


viscous  solution  is,  in  this  case,  4-  =  consteint  every^rtiere,  i.e.,  the 


frictionless  notion  of  the  walls  does  not  affect  the  fluid  In 


ence  value  on  the  obstacle  and  on  a  horizontal  line  through  the  center 


of  the  obstacle*  Along  some  line  Joining  the  vails,  but  not  Including 


As  a  final  example  of  an  Initial  solution  ve  simply  attach 


nocUfying  fit  dowtwaPd  should  conditions  arise  In  which  these  criteria  “i.J  calculation  Is  discussed  in  detail  In  Section  6-C.  The  ♦  calcul 
may  be  violated.  There  is  no  need  for  upward  revision  of  6t  since  Its  advanced  cu  values,  is  discussed  in  Section  6-D.  Ve- 
slze  does  not  strongly  Influence  the  speed  of  calculation  (see  Section  J).  locity  components  corre^ndlng  to  new  ♦  values  arc  calculated,  and 


Set  Index 


for  each  mesh  point  to  determine  the  treatment  to  he  given  the  vorticity 


properties  are  satisfied, 


A  CD 


except  that  u  value  Just  Inside  the  vallc 


That  is 


We  begin  vlth  considerablonB  of  the  faaniliar  Bynmetrlc  vake  for  order  of  Reynolds  numbers  given  in  Fig,  7.  In  all  cases  the  integrated 

low  Reynolds  numbers  (see  Fig.  7).  The  Reynolds  number  for  the  flows  vorticity  magnitude  grew  from  that  initially  present,  having  been  the 


obstacle.  For  exan^jle,  the  integrals  of  the  magnitude  of  vorticity  in  the  ir.ally.  We  assume  that  once  circular  flow  has  become  established,  vor 

field  at  the  specified  time  are  in  ratios  of  about  20:23:26:36,  in  the  ticity  tends  to  become  trapped  into  regions  of  extreme  values  near  the 


ation,  since  less  tran^ort  of  vorticlty  can  occur  We  proceed  to  a  consideration  of  a  flow  sequence  (Figs*  8  and  9)  of 

such  points.  This  would  happen  more  readily,  of  course,  if  the  a  vortex  street  formation  for  the  case  R  =  100.  The  numerical  program  ha 

diffusion  rate  were  slow.  In  all  cases  observed  the  detached  extremes  been  shown  to  persist  in  giving  a  synmetrlc  flow  pattern  unless  a  pertur- 

ortlcity  diffused  away,  forming  less  peaked  extremes.  Indications  bation  is  introduced.  We  introduce  a  perturbation  by  augaentlng  the  vor- 
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tine  R  has  dropped  to  ab<mt  hali'  Its  original  value. 
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be  augmented  by  the  mainstream.  Continuation  of  this  process  leads  discussion  we  have  said  nothing  of  the  vorticity  and  Its 

fluence  on  what  occurs.  It  is,  of  course,  the  mathematical  basis  for 


formation  of  the  eddies,  hut  the  motion  of  the  fluid 


owing  to  diffusion,  vdiile  the  integrated  vorticity  associated  with  the 


eddy  increases  through  the  addition  of  vorticity  tran^rted  from  the  ob- 


average  flow  rate  has  gone  down  and  the  new  eddy  cannot  again 


Since  the  velocity  coDi)onent  values  are  given  everywhere  in  the 


centers  of  i];'  contours  are  at  a  greater  distance  from  the  street  center  again  note  that  the  tJt'  extreme  is  displaced  outward  from  the 

line  than  are  extremes  in  vorticity,  in  agreement  with  the  discussion  fixed  wall  (or  the  center  line)  relative  to  the  vorticity  extreme, 

given  hy  Shaefer  and  Eskinazi,^^  Because  of  the  weakness  of  the  eddies  after  the  first,  only  the  first 

is  well  portrayed  on  the  plot. 


components  are  of  opposite  sign,  then  we  would  have  to  let  = 


ing  diagram  and  also  that  added  region  which  has  been  found  to  exist 


con^Jutation  time.  It  is  for  this  reason  that  other  possible  difference  ^  first  exainple  of  the  right-hand  column  we  include  a  case 

forms  have  been  sought. 


,  according  to  the  analysis,  in  a  stable  region.  We  account  that  the  v  =  0  case  is  conjjletely  unstable.  If  this  case  could  be  made 

growth  in  that  the  perturbation  is  not  sufficiently  small  for  stable,  then  the  stability  properties  would  be  lugaroved  for  all  values 

analy  to  be  valid  in  this  near  proximity  to  the  unstable  region.  of  v.  We  therefore  look  for  a  difference  form  of  the  vortlclty  equation 
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and  Harlow,  F.  H.,  Phys*  Fluids,  to  te  issued, 
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Fig,  1  •  Layout  of  the  lattice  of  points,  indicating  where  variables  i 
defined  relative  to  the  mesh  and  the  geometry  of  the  channel, 


r- 

o 


60 

£ 


I 

VO 

O 


Fig.  10.  First  part  of  streakllne  sequence  for  R 


Qualitative  comparison  between  exainples  of  e^^perimental  and 
calculated  streakline  patterns.  Upper  frame  obtained  by  Thom, 


Fig.  17.  Second  part  of  streaJcUne  sequenci 


Fig,  19*  Second  part  of  streamline  sequence 


Fig,  21 •  Transforaed  streamline  sequence  for 


Fig.  22.  First  part  of  streamline  sequence  for  R  =  1000. 

d/h  =  l/6.  Perturbation  introduced  at  t  ~  0,14 


:nB.  24.  streakline  sequence  for  R  =  1000.  streamline  sequences  for  R  1000  and  R 


turbatlon  introduced 


Fig.  32.  First  part  of  streakline  sequence  for 


Fis-  57«  Transformed  streamline  sequence  for 


Fig.  58*  First  part  of  streakline  sequence  for  R  =  500.  Fig.  59.  Second  part  of  streakline  sequence  for 


Fig.  1*2.  First  part  of  streakllne  sequence  for  R  =  1(XX3.  Fig.  4-5.  Second  part  of  streakline  sequence  for 


projection 


Fig*  47.  Second  part  of  the  streakline  sequence  of  flow  about  a  pro¬ 
jection  from  a  fixed  wall. 


Various  representations  of  the  flow  past  an  obstacle  for  Pie.  49.  Various  representations  of  the  flow  about  a  projection 

=  ’/la,  Au=  6.4,AP=  1/3,  .  fixed  wall.  At  =  1 /6,  At'  =  l/6,  Au)  ^eft)  =  6.4, 

d/h  -  1/u.  Perturbation  introduced  at  t  ~  0.20.  Am  (’■rleht')  =  3.2.  AP  =  1.6. 


Fig.  50.  Patterns  of  jmstable  growths  of  a  perturbation  using  a  nontime- 


